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Abstract 

We establish intra-individual and inter-annual variability in European badger (Meles meles) autumnal nightly activity in 
relation to fine-scale climatic variables, using tri-axial accelerometry. This contributes further to understanding of causality in 
the established interaction between weather conditions and population dynamics in this species. Modelling found that 
measures of daylight, rain/humidity, and soil temperature were the most supported predictors of ACTIVITY, in both years 
studied. In 2010, the drier year, the most supported model included the SOLAR*RH interaction, RAIN, and30cmTEMP 
(w = 0.557), while in 2012, a wetter year, the most supported model included the SOLAR*RH interaction, and the 
RAIN*10cmTEMP (w = 0.999). ACTIVITY also differed significantly between individuals. In the 2012 autumn study period, 
badgers with the longest per noctem activity subsequently exhibited higher Body Condition Indices (BCI) when recaptured. 
In contrast, under drier 2010 conditions, badgers in good BCI engaged in less per noctem activity, while badgers with poor 
BCI were the most active. When compared on the same calendar dates, to control for night length, duration of mean badger 
nightly activity was longer (9.5 hrs ±3.3 SE) in 2010 than in 2012 (8.3 hrs ±1.9 SE). In the wetter year, increasing nightly 
activity was associated with net-positive energetic gains (from BCI), likely due to better foraging conditions. In a drier year, 
with greater potential for net-negative energy returns, individual nutritional state proved crucial in modifying activity 
regimes; thus we emphasise how a 'one size fits all' approach should not be applied to ecological responses. 
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Introduction 

In seasonal ecosystems, biological and chemical processes [1], 
[2], [3], along with phenological responses (e.g., see: [4]), generally 
respond to annual cyclicity in primary productivity [5], weather 
patterns, and hydrological, nutrient, and carbon cycles [6], [7], 
[8] . A species' climatic niche, and range of tolerance around this 
niche, are thus defining elements in evolution [9], [10], [11]. 
Considering that ecological responses are typically non-linear [12], 
[13], even minor changes in environmental conditions can have 
substantial effects across trophic levels [14]. For consumers, 
seasonal variation in food supply and associated thermal stresses 
while foraging create energetic demands [15], [16], [17], driving 
physiological and/or behavioural adaptations - for example, 
endocrinological changes [18], fat deposition [19], [20], [21], 
winter hibernation/torpor [22], [23], life-cycle synchronicity [24], 
[25], migration [26], and food caching [27], [28]. Matching 
species' mechanistic responses to fine-scale environmental vari- 
ability, while imperative, proves challenging (e.g., [29], [30], [31]). 

To these ends, recent advancements in accelerometery 
technology capable of measuring both animal orientation and 
the dynamics of movement, have considerably enhanced the 
capacity to achieve behavioural activity sampling on a fine-scale 
[32], [33], [34], [35]. Accelerometry has elaborated a wide range 



of behaviours - such as diving [36], [37], feeding [38], [39], and 
mating behaviour [40] - which often cannot be recorded as 
systematically by other techniques in wild, cryptic animals [41], 
[42]. Important in this study, the activity data accelerometers 
generate can be utilised to provide an accurate proxy for energy 
expenditure [43]. 

Here, we apply accelerometry as a tool to investigate the 
interaction between proximate weather conditions and European 
badger (Meles meles) activity regimes ([44], [13]); that is, this 
approach enables us to assess how activity is promoted or inhibited 
by micro-climatic conditions 'per individual'. 

The European badger provides a particularly good model for 
studying the effects of weather, because of its sensitivity to climatic 
conditions ([45], [20], [44], [13]) and its trophic specialism, 
favouring earthworms, Lumbrkus terrestris, in lowland Britain and 
Ireland [46], [47], [48], [49]. 

Badgers undergo autumnal weight gain ([50], [20]) in advance 
of a flexible extent of winter torpor [51], where they rest inside 
extensive dens, termed setts [52], [53], [54]. This inactivity period 
is most evident in higher latitude regions of their range, where 
winter frosts restrict earthworm availability ([55], [56], [57], [58]). 
Earthworms become less available under frosty conditions [59], 
and so in order to maintain or increasing their body-condition in 
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the lead up to food supply restriction (i.e., over the autumn) 
badgers must attempt to balance trade-offs between the energetic 
expenditure of their activity (foraging and non-foraging), against 
energy conserved by inactive periods within their setts. As a 
consequence early autumn - specifically September - weather 
conditions present an important transition point in badger activity 
regimes influenced by variability in the change over from 
'productive' to 'less/non-productive' conditions [20], [44], where 
more activity has been observed under sub-optimal conditions, to 
attempt to compensate for lower foraging success [60], [58], [61]. 

In addition to nutrition, there are also other individual-specific 
motivators of activity, such as territorial marking [62], [63], scent 
communication [64], and associated social contact networks [65] 
potentially pertinent to reproductive success and thus fitness [66] . 
Consequendy, complex individual-specific differences in activity 
patterns are apparent [67], [68], [69]. 

Given that interactions between badger population dynamics 
and weather metrics have been established [20], [44], [13], here 
we aim to advance understanding causality by examining how 
variation in autumnal weather conditions in two non-consecutive 
years influenced activity budgets. We use overall dynamic body 
acceleration (ODBA) from tri-axial accelerometers [41], [70], to 
expose real-time variation in badger activity levels through late 
summer/early autumn. We also postulate that variation in 
nutritional state will yield differing individual responses to 
environmental stressors, evidenced as differences in activity levels. 

Methods 

The study was conducted on a high-density badger population 
(36.4 badgers/km 2 , SE= 2.55 badgers/km 2 [71]) in Wytham 
Woods, Oxfordshire (GPS reference: 51°46'26"N; 1°19'19"W); a 
424-ha site of mixed semi-natural woodland (for detailed 
description of the study site see [20], [72]). As part of a trapping 
regime undertaken since 1987, badgers were caught in box traps 
(85x37x38 cm), baited with 150 g of peanuts (but without pre- 
baiting, making the overall effect of food supplementation 
inconsequential - see [71]) placed near the entrances of setts 
and outlying holes (outliers) over two weeks in June (spring), 
September (summer) and November (autumn), with one week of 
trapping in January (winter) in some years [71]. Captured animals 
were transferred to holding cages, between 06:30h and 08:00h the 
following morning, and transported to a field station using an 
ATV and trailer. At the field station, individuals were sedated by 
an intramuscular injection of ketamine hydrochloride (0.2 mL/kg 
body weight) [73], [74], sexed, measured to the nearest mm (tip of 
snout to base of sacrum, laid dorsally), and weighed to the nearest 
100 g. 

From these morphometries, we used the ratio of body length to 
weight (L:W) to calculate Body condition indices (BCI) (log mass 
vs. log length), as a measure of an individual's nutritional state 
[75], [76]. To monitor seasonal and individual-specific changes in 
BCI, we calculated the index for the summers in which animals 
were instrumented and, provided the animals were (re-)captured, 
the spring and autumn preceding/following these summers. 

We instrumented four badgers in 2010 (3 males and 1 female) 
and four different individuals in 2012 (3 males and 1 female) with 
prototype tracking tags [39 mm (l)x22mm (w)xl2mm (h)] 
mounted on a leather collar, with a complete assembly weight of 
105 g (for detailed description of the collars and data collection 
system see: [41], [70]). Henceforth individual badgers were 
identified by their unique collar ID. Tags were equipped with 
tri-axial accelerometers, sampling pitch, yaw, roll at 8 Hz, as well 
as radio-frequency identification (RFID), and magneto-inductive 



sensors (MI) - where a receiver measured the strength of magnetic 
fields generated by multiple, low frequency loop antennas, 
allowing the position of the collared animal to be determined in 
3-D (see: [41], [70]). 

After handling and collaring, badgers were held in a recovery 
room for 3 h, before being released at the locations where they 
were captured. Data were compressed and stored onboard the 
collars, and communicated wirelessly over a 2.4 GHz 802.15.14 
link when the animal came above ground, making recapture for 
data collection unnecessary. In 2010, data were collected from 
September 28 until batteries ran out on November 24, i.e., for a 
total of 58 days. All batteries lasted for the same duration except 
for those in one collar (N30), which ran out after only 45 days. In 
2012, data were collected from September 1 until all batteries ran 
out on October 30, i.e., for a total of 60 days. 

We converted these acceleration values for the three axes/ 
channels to overall dynamic body acceleration (ODBA), using the 
method developed by Wilson et al. [77] . That is, each channel was 
smoothed individually; using running means over 10 min time 
frames. The values for these smoothed data, for any particular 
time frame, were then subtracted from the corresponding 
unsmoothed data for that time frame. This gives a value for g 
resulting from dynamic acceleration that we converted into 
absolute positive units using the LI norm, combining values from 
all three channels, to produce an ODBA metric. ODBA values 
were compared to a threshold, determining badger activity from 
inactivity for each data period (ACTIVITY), (see Appendix SI, 
Figure SI, Table SI). To satisfy assumptions of normality, hourly 
measures of activity were subjected to angular transformations for 
all statistical analyses. 

Meteorological records were obtained from the UK Environ- 
mental Change Network (ECN)'s weather station in Wytham 
Woods (T08) and are availably freely. In order to separate the 
proximate effects of weather on badger ACTIVITY from inter- 
annual variation, we decomposed climatic data into 7 environ- 
mental variables, known to influence badger behaviour and 
physiology [44], [45], [52], [54]: (1) solar radiation, providing a 
measure of photoperiod (SOLAR, in W/ m 2 ) with the additional 
benefit of including cloud cover and thus more subtie influences on 
the timing of dusk emergence; (2) relative humidity (RH, in %); (3) 
mean air temperature (TEMP, in °C); (4) total rainfall (RAIN, in 
mm); (5) soil temperature at a depth of 10 cm (lOcmTEMP, in °C); 
(6) soil temperature at a depth of 30 cm (30cmTEMP, in °C); (7) 
volumetric water content of the soil (WATER, mVm 3 ). We then 
modelled the effects of these climatic components on badger 
ACTIVITY both within each study period and between years.. 
Given the importance of photoperiod on badger activity [52], we 
utilised the terms 'per diem' and 'per nocterri to distinguish between 
daylight and night-time hours respectively, and avoid confusion 
associated with the dual meaning of the term 'per day'. 

Statistical Analyses 

We performed a preliminary analysis of daily measurements 
using analyses of variance (ANOVA) with analyses of covariance 
(ANCOVA), paired-sample t-tests, coefficients of variation (CV), 
and Levene's tests (p = 0.05). As ACTIVITY data were collected 
over two discrete windows within the two years, analyses relating 
climate to ACTIVITY were performed for each period separately 
and then compared. 

We then applied an information theoretic approach [78], [79] 
to examine finer scale-hourly data, using Akaike's information 
criterion (AIC), in a mixed model approach. SOLAR was included 
in all multi-variate models to reflect the nocturnality of badgers, 
and to avoid superfluous testing; no model included more than one 



PLOS ONE | www.plosone.org 



2 



January 2014 | Volume 9 | Issue 1 | e83156 



Autumnal Activity of the European Badger 



measure of temperature (i.e., TEMP, lOcmTEMP, or 
30cmTEMP). In addition, because badgers differed individually 
in their activity regimes, we included individuality as a random 
effect variable. We acknowledge that although a degree of 
variation in activity regimes was likely due to inter-sexual 
variation, sample sizes here did not permit us to investigate this 
in detail. Given that seven climatic variables were derived for each 
year, we constructed 26 models without interaction and 33 models 
that included (an) interaction term(s). We therefore evaluated 59 
models with climatic parameters as predictors of badger ACTIV- 
ITY. For each model, we derived the AIC, which we used to rank 
the support for each model (a lower value indicating greater 
support), as well as the delta AIC (A;) in relation to the highest 
ranking model and the Akaike weight (w) [78]. As per Anderson 
[80], A; cut-off points were not used, though models with A ; values 
above 9-1 1 were considered to have relatively little support, and 
models with A; values >20 had no empirical support [79]. 
Statistical analyses were performed using RStudio v. 0.96.331. 
Although statistical analyses were performed on transformed data, 
for visual purposes, we used the un- transformed data in all figures. 

Ethics 

All trapping and handling procedures were subject to ethical 
review and were performed under Badger Act (1992) licence 
(currently 20104655) from Natural England and UK Animals 
(Scientific Procedures) Act, 1986 licence from the Home Office 
(currently PPL30/2835). 

Results 

Autumnal Climatic Conditions 

Between the autumn study periods in the two years, we found 
no difference in daily TEMP (ANOVA: F [M1B] = 1.692, 
p = 0.196), (summarised in table 1), nor between TEMP variability 
(2010: CV 27.8%, n = 58; 2012: CV = 27.1%, n = 60; Levene's: 
F[i,i i6] = 0.662, p = 0.432). Badgers are nocturnal however, and, 
controlling for night length, 2012 proved significantly colder (daily 



TEMP 4°C cooler) than 2010 (paired sample t-test: t, 2 = -8.352, 
p<0.0005) over matched calendar dates. 

Decomposing temperature variables further, we found that in 
2012 TEMP decreased steadily over the duration of the study 
period. In contrast, in 2010 there was a cold interlude of 16 days in 
early November (Nov. 7-Nov. 22), where TEMP decreased 
rapidly, from 16°C per day to 8.5°C over 5 days, and averaged 
7.2°C±1.8 SD over the length of the interlude compared to 
13.4°C± 1.8 SD over the first 40 days of the study period. TEMP 
variability over this interlude, through to the end of the study 
period, was also significantly more variable than the first month 
and a half of the study (CV = 29.5%, n= 18; CV= 14.9%, n = 40 
respectively; Levene's: F[ lj56 ] =8.614, p = 0.005) (Fig. 1). 

As with TEMP, there was no evidence for a difference in daily 
lOcmTEMP (ANOVA: F [1>116] = 0.384, p = 0.536), but 2012 was 
significandy colder than 2010 when compared directly on the 
same calendar dates (paired sample t-test: t 57 = 13.546, p<0.0005) 
with mean daily lOcmTEMP being 3.5°C cooler. Similarly, there 
was no evidence for a difference in 30cmTEMP, (ANOVA: 
F[i,ii6] = 0.003, p = 0.995), but 2012 was significantly colder than 
2010 when compared on matched calendar dates (paired sample t- 
test: t 57 = 24.707, p<0.0005) with mean per day 30cmTEMP 
being 2.6°C colder. Within years, lOcmTEMP was significandy 
more variable than 30cmTEMP (2010: CV = 49.36%, n = 58; 
CV= 20.21%, n = 58 respectively; Levene's: F [1>114] = 44.499, p< 
0.0005; 2012: CV= 18.77%, n = 60; CV= 13.41%, n = 60 
respectively; Levene's: Fp ns] =4.781, p = 0.031). Between years, 
both lOcmTEMP and 30cmTEMP in 2010 were significandy 
more variable than in 2012 (Levene's: lOcmTEMP; 
F[i, n6] =51.119, p<0.0005; 30cmTEMP: F [U16] = 7.01 1, 
p = 0.009). 

With respect to measures of humidity, although there was no 
evidence for a difference in RAIN between these years (ANOVA: 
F [W i 6] = 1.089, p = 0.299), the 2012 study period had significandy 
higher RH (ANOVA: Fp.ne] = 7.286, p = 0.008), and WATER 
(ANOVA: ^[1,116] = 20.841, p<0.0005) than 2010. 



Table 1. Climatic parameters and ACTIVITY. 





2010 


2012 


TEMP (°C) 


11.75±3.24 


11. 00 ±2.96 


TEMP paired TO 


13.34±2.11 


9.41 ±2.24 


lOcmTEMP TO 


13.07 ±2.62 


12.78±2.38 


lOcmTEMP paired (=C) 


14.65 ±1.44 


11.11 ±1.38 


30cmTEMP (°C) 


13.91 ±1.79 


13.90±1.84 


30cmTEMP paired (=C) 


15.05 ±0.86 


12.47±0.82 


RH (%) 


83.80±7.43 


87.84±8.64 


RH paired (%) 


83.09±8.48 


9 1.67 ±6.45 


RAIN (mm) 


1.82±3.35 


2.73±5.75 


RAIN paired (mm) 


2.21 ±3.81 


2.98±5.68 


WATER (m3/m3) 


0.306±0.118 


0.383±0.053 


WATER paired (m3/m3) 


0.301 ±0.1 17 


0.42±0.019 


ACTIVITY (hrs) 


9.48±3.32 


8.27 ±1.86 


ACTIVITY paired (hrs) 


11.25±1.02 


9.22±1.17 



Mean values for climate parameters and ACTIVITY ± SD over the duration of the 
study periods in 2010 (58 days) and 2012 (60 days). Paired parameters represent 
means over the same calendar dates (Sep. 28-Oct. 30). 
doi:10.1371/journal.pone.0083156.t001 
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20 1 J 



f «,<> V * 



Aug/30 Sep/13 Sep/27 Oct/11 Oct/25 Nov/8 Nov/22 

Figure 1 . Mean daily TEMP over the study periods. Mean daily air 
temperature (TEMP) over the course of the study periods in 2010 and 
2012. Between the two years, we found no difference in daily TEMP 
(ANOVA: Fn,,, 6] = 1.692, p = 0.196), nor between TEMP variability (2010: 
CV 27.8%, n = 58; 2012: CV = 27.1%, n = 60; Levene's: F [1-116] = 0.662, 
p = 0.432) but 2012 proved significantly colder (daily TEMP 4°C cooler) 
than 2010 (paired sample t-test: t 32 = -8.352, p<0.0005) over matched 
calendar dates. 

doi:10.1371/journal.pone.0083156.g001 
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Badger Activity 

Inter-annual differences in badger activity. From Sep- 
tember through October in 2012 per noctem ACTIVITY increased 
with night length and thus correlated negatively with photoperiod 
(SOLAR: (ANCOVA: F [lj230] = 74.31 1, p<0.0005). In 2010, 
however, ACTIVITY was monitored over a later period, from 
October through November, and per noctem ACTIVITY decreased 
with lengthening nights, resulting in a positive correlation with 
SOLAR (ANCOVA: F [1216] = 16.955, p<0.0005). Over the study 
period in both years, there was no significant difference between 
individuals (ANCOVA 2010: F [li216] = 3.301, p = 0.071; 2012: 
F[i,230] = 3.578, p = 0.060) This suggested that some other 
variable, such as weather, was counteracting night length. Badgers 
were significantly more active in 2010 than in 2012 (ANOVA: 
F[i,i454] = 13.510, p<0.0005) (Table 1), when compared on the 
same calendar dates between years, to control for night length. As 
well as being more active, the duration per noctem ACTIVITY was 
significandy more variable in 2010 than it was in 2012 
(CV = 35.3%, n = 58; CV = 22.7%, n = 60 respectively; Levene's: 

F[l,116] ~ " 



3.462, p = 0.004) (Fig. 2). 



Interaction between activity and climatic variables. To 

explore how the weather conditions influenced per noctem 
ACTIVITY between years specifically, we determined the 
relationship between the climatic variables and ACTIVITY using 
ANCOVAs. The study period in 2010 was comparatively drier 
than 2012, with a pronounced cold interlude, and more variable 
soil temperatures. In 2010, we observed no correlation between 
mean badger ACTIVITY and RH (Fig. 3 a) (F [1>216] = 0.458, 
p = 0.50). There also was no correlation between activity and 
RAIN (F [lj2 i 6 ] = 0.088, p = 0.77), but a positive correlation 
between ACTIVITY and TEMP (Fig. 4 a), lOcmTEMP, 
30cmTEMP, and WATER (F [1>216] = 86.493, p<0.0005; 
F[i j2 i6] =8.942, p<0.005; F [1>216] = 120.510, p<0.0005; 
F[i,2i6] = 7.018, p = 0.008 respectively). In all cases, there was no 
significant difference between individuals (p&0.063). 

As well as being wetter than 2010, soil temperatures were also 
less variable during the 2012 study period, and we found that 
ACTIVITY correlated significantly and positively with RH 
^[1,230] = 102.593, p<0.0005) (Fig. 3 b), WATER 
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Figure 3. Mean per noctem ACTIVITY and RH over the study 
periods. Mean per noctem ACTIVITY and relative humidity (RH) over 
the course of the study periods in 2010 (A) and 2012 (B). In 2010, there 
was no correlation between ACTIVITY and RH (F M , 216] = 0.458, p = 0.50). 
In 2012, ACTIVITY correlated significantly and postively with RH 
(F n ,230] = 1 02.593, p<0.0005). 
doi:10.1371/journal.pone.0083156.g003 
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Figure 2. Mean per noctem ACTIVITY over the study periods. 

Mean per noctem ACTIVITY over the course of the study periods in 2010 
and 2012. The duration of per noctem ACTIVITY was significantly more 
variable in 2010 than it was in 2012 (CV = 35.3%, n = 58; CV = 22.7%, 
n = 60 respectively; Levene's: F [q>q , 6] = 8.462, p = 0.004), this was also true 
when compared on the same calendar dates, to control for night 
length, significantly longer (ANOVA: ^,454]= 13.510, p<0.0005) in 
2010. 

doi:1 0.1 371 /journal.pone.00831 56.g002 



(F[i j2 30] = 71.524, p<0.0005), and RAIN (F [1)230] = 13.368, p< 
0.0005). In addition, we found a negative correlation between 
ACTIVITY and lOcmTEMP, and 30cmTEMP (F [1;230] = 22.205, 
p<0.0005; F [li230] =62.265, p<0.0005 respectively), but no 
correlation between mean ACTIVITY and TEMP (Fig. 4 b) 
(F[i ; 230] = 2.456, p = 0.1 18). Again, in all cases, there was no 
significant difference between individuals (p&0.062). 

Modelling Activity-climate variable interactions . Informed 
by our observation that the relationship between these environ- 
mental variables and ACTIVITY varied between years, we used 
information-theoretic models, and finer-scale hourly measures, to 
ascertain the strongest predictors of per noctem ACTIVITY. For 
2010, the most supported model included the SOLAR*RH inter- 
action, RAIN, and 30cmTEMP (w = 0.557), however a model that 
included the SOLAR*RH interaction, and the RAIN*30cmTEMP 
interaction also had some support (w = 0.442). A; values indicated 
too much information loss to warrant consideration of other models. 
For 2012, the most supported model included the SOLAR*RH 
interaction, and the RAIN*10cmTEMP interaction (w = 0.999). 
Again, A; values indicated too much information loss (Table 2) to 
warrant consideration of other models. Thus, while our uni-variate 
analyses revealed that the relationship between each environmental 
variable and ACTIVITY varied inter-annually, this modelling 
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Figure 4. Mean per noctem ACTIVITY and TEMP over the study 
periods. Mean per noctem ACTIVITY and daily air temperature (TEMP) 
over the course of the study periods in 2010 (A) and 2012 (B). In 2010 
ACTIVITY correlated significantly and postively with TEMP 
(F[i,2i6] = 86.493, p<0.0005), where in 2012, there was no correlation 
between ACTIVITY and TEMP (F [1 , 230 ] = 2.456, p = 0.1 18). 
doi:10.1371/journal.pone.0083156.g004 

showed that measures of daylight, rain/humidity, and soil 
temperature were the most supported predictors of ACTIVITY in 
both years. 

Individual Differences in Badgers Activity Regimes 

Within this climatic framework, in both years ACTIVITY also 
differed significandy between individuals (ANOVA: 
F [3j5 280] = 4.866, p = 0.002; F Pi 5 496] = 3.382, p = 0.017 respective- 
ly). We therefore compared BCI, as a proxy of nutritional state, 
against individual activity regimes, as well as the population mean 
(PM) values for BCI ± the standard deviation (SD). 

In 2010, badgers with higher BCI, exhibited less per noctem 
activity than badgers with lower BCI. When first captured in 
September, only one badger (N30) had a significandy poorer BCI 
than the PM. This individual also exhibited the longest mean 
ACTIVITY duration per noctem compared to all other badgers 
tracked that year, being active on average 30 min longer than 
conspecifics; this greater activity was statistically significant when 
compared to badgers N3 1 and N34 (Tukey post-hoc: p = 0.003; 
p = 0.020 respectively). N31 and N34 had BCIs significandy 
higher than the PM. N33 was, on average, also active for 30 min 
longer per noctem than N34, however this was not statistically 



significant (ANOVA: F [lj27!!9] = 2.801, p = 0.094). When recap- 
tured in November, N33 had a lower BCI than N34, and was 
within the PM, while N34's BCI had increased, and was still well 
above the PM (Fig. 5a). 

In contrast in 2012, badgers with the longest per noctem activity 
exhibited higher BCIs when recaptured in November. BCI for 
N48 and N49 both exceeded the PM when collared in September. 
When recaptured in November, only N48 remained above the 
PM, and had been more active (on average 20 minutes longer per 
noctem) than N49, for which BCI had fallen to within the PM. 
Badger N44 had a lower BCI than both N48 and N49, and was 
below the PM in both September and November. On a per noctem 
basis, N44 was, on average, active for 40 min less than N48 and 
20 min less than N49. Badger N45 was the least active over the 
course of the study period, and was within the PM when collared 
in September, however it was not recaptured in November, and so 
no conclusions can be drawn for that individual concerning 
autumnal changes in BCI (Fig. 5b). 

Discussion 

Badgers are responsive to climatic conditions [45], [20], [44], 
[13] and, in this study population, autumnal conditions have been 
found to present a critical period for all aspects of badger 
population dynamics [44] , testing the ability of badgers to balance 
energy budgets. Applying fine-scaled measures of activity allowed 
us to determine with high resolution that individual badger activity 
regimes varied in response to proximate weather, and with respect 
to body-condition. These clear inter-annual differences in the way 
individuals responded to varying climatic conditions evidence that 
biological responses, phenology and cyclicities in annual physio- 
logical and behavioural regimes should not be assumed to follow 
the same course each year. 

Although there are various factors that can influence badger 
activity (for example Goransson [81] found that, in autumn, 
badgers spent over 2-5 h daily digging and collecting material for 
bedding), our modelling indicates that interactions between solar 
radiation and humidity, rainfall, and soil temperatures were 
significant drivers of badger per noctem autumnal activity. 

In relation to Nouvellet et al's [13] observation that survival 
rates decline outside of the typical annual range of daily predicted 
temperature norms in this population, we found that badgers were 
significandy more active in the drier year (2010) than a wetter one 
(2012). 

Furthermore, we identified the generality for those badgers 
attaining the highest BMI to also be the most active individuals in 
2012; the wetter of the two years. This implies that whatever 
specific activities contributed to their nightly regime, foraging must 
have being among these, and returned net-positive energetic gains. 
In 2010, with drier conditions, we observed a more subtle 
response: 'thinner' badgers (BCKPM) were more active than 
'fatter' badgers (BCI>PM). 

Given that badgers in this region [46], [47], [48], [49], [82], as 
well as in others (e.g., [83], [58]), favour earthworm prey (see 
above), and given that earthworms surface and become most 
available to predators under mild, moist microclimatic conditions 
[57], [84], [58], this syllogism leads us to strong inference that 
activity relates in major part to foraging success. Badgers with 
higher BCI prove freer to obviate the risk of net-negative energy 
returns under less optimal conditions by reducing their nighdy 
activity levels. 

Our tri-axial accelerometry results resonate with other studies 
that report increased badger foraging-related activity under sub- 
optimal climatic conditions (e.g., in Spain [60], Poland [58], and 
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Figure 5. Variation in individual nutritional states from the population. Individual measures of nutritional state (L:W) in comparison to the 
population mean (PM) ± the standard deviation (SD) in June, September, and November in 2010 (A) and 2012 (B). Red points indicate badgers with 
comparatively long ACTIVITY, where blue points indicate badgers with comparatively short ACTIVITY. 
doi:1 0.1 371 /journal.pone.00831 56.g005 



Switzerland [61]). Under conditions that limit earthworm 
availability (i.e., below 2°C: [55], [56], [57], [58]), Kowalczyk et 
al. [58] report that badgers attempt to compensate by devoting 
over 2.5 hrs more to foraging. In addition, maintaining body- 
temperature is more energetically expensive when active in cold, 
wet conditions [85], [86] - see also [44], [13]. This defines a suite 
of optimal weather conditions that can influence badger foraging 
success [87], ability to improve /main tain body condition [20] and 
thus ultimately survival probability [44], [13]. 

Based on the inter-individual differences in activity regimes we 
observed, we emphasise that a 'one size fits all' approach cannot 
be applied to ecological responses. State-dependent risk-taking is 
well documented across a diversity of taxa [88], where changes in 
behaviour to maintain positive energy budgets have been linked to 
climate variability (e.g., caribou (Rangifer tarandus) [89]; badgers 
[44]; Eurasian beavers (Castor fiber) [4]). According to Clark [90], 
individuals with greater reproductive value - a function of body 
condition in badgers [91], [20], should increase risk aversion in 
comparison to individuals with less reproductive value, in order 
not to jeopardise survival and future reproduction. Conversely, 
some individuals will take greater risks to achieve future fitness 
returns than individuals with greater reproductive value [92] . We 
infer that predictable conditions seem to be important for optimal 
survival dynamics, allowing individual badgers to prospect risk 
most effectively [93] and optimise their annual routines [94]. 

In conclusion, these accelerometry data further advance 
understanding of badger-climate interactions, with detail allowing 
us to begin to dissect the causal mechanisms involved. This level of 
resolution is essential - beyond identifying indicative trend 
relationships - in order to identify the ecological consequences 
of weather patterns (or indeed other influential factors in other 



systems) and understand individual-specific adaptive responses, in 
turn better informing ability to devise conservation policy [95] . 

Supporting Information 

Figure SI Perceived activity as a function of the 
threshold value. Percent active as a function of the threshold 
value (T) for badger N44 over the course of the study period. The 
threshold value chosen for this individual was 10. 
(TIF) 

Table SI Perceived hourly activity. Hourly activity (in 
minutes) for four collared badgers over the study period in (A) 
2010 and (B) 2012. For queries on finer-scale data please contact 
the corresponding author. 
(XLSX) 

Appendix SI 

(DOCX) 

Acknowledgments 

This study would not have been possible without the efforts of the badger 
trapping team, past and present. We thank the UK Environmental Change 
Network (ECN) for the provision of meteorological data. We are grateful to 
Paul Johnson for instruction in analytical procedures. 

Author Contributions 

Conceived and designed the experiments: AM NT CN. Performed the 
experiments: CDB AM SAE. Analyzed the data: MJN. Contributed 
reagents/materials/analysis tools: NT AM. Wrote the paper: MJN CN 
DWM. 



References 



1. Gaedkc U, Strailc D (1994) Seasonal changes ol the quantitative importance of 
protozoans in a large lake. An ecosystem approach using mass-balanced carbon 
How diagrams. Marine Microbial Food Webs 8: 163-188. 



2. Catalan J, Ventura M, Camarero L, Granados I, Thics H, ct al. (2002) Seasonal 
ecosystem variability in remote mountain lakes: implications for detecting 
climatic signals in sediment records. J Paleolimnol 28: 125-46. 



PLOS ONE | www.plosone.org 



7 



January 2014 | Volume 9 | Issue 1 | e83156 



Autumnal Activity of the European Badger 



3. Ishii HT, Tanabc S, Hiura T (2004) Exploring the relationships among canopy 
structure, stand productivity, and biodiversity of temperate forest ecosystems. 
Forest Science 50: 342-355. 

4. Campbell RD, Newman C, Macdonald DW, Rosell F (2013) Proximate weather 
patterns and spring green-up phenology effect Eurasian beaver (Castor fiber) 
body mass and reproductive success: the implications of climate change and 
topography. Global Change Biology 19: 1311-1324. 

5. Nemani RR, Keeling CD, Hashimoto H, Jolly WM, Piper SC, et al. (2003) 
Climate-driven increases in global terrestrial net primary production from 1982 
to 1999. Science 300: 1560-1563. 

6. BarthJAC, VcizerJ (1999) Carbon cycle in the St. Lawrence aquatic ecosystems 
at Cornwall (Ontario), Canada: seasonal and spatial variations. Chem Gcol 159: 
107-128. 

7. Huxman TE, Turnipseed AA, Sparks JP, Harley PC, Monson RK (2003) 
Temperature as a control over ecosystem C02 fluxes in a high-elevation, 
subalpine forest. Oecologia 134: 537-546. 

8. Greco S, Baldocchi DD (2006) Seasonal variations of C02 and water vapour 
exchange rates over a temperate deciduous forest. Global Change Biology 2: 
183-197. 

9. Caswell H (1982) Optimal Life Histories and the Maximization of Reproductive 
Value: A General Theorem for Complex Life Cycles. Ecology 63: 1218-1222. 

10. Thuiller W, Lavorel S, Araulo MB (2005) Geographical extent and niche 
properties as predictors of species sensitivity to climate change. Global Ecology 
and Biogeography 14: 347-357. 

11. Wolf M, van Doom S, Lcimar O, Weissing FJ (2007) Life-history trade-offs 
favour the evolution of animal personalities. Nature 447: 581-584. 

12. May RM (1986) When two and two do not make four: nonlinear phenomena in 
ecology. Proc R Soc Lond B 228: 241-266. 

13. Nouvellet PM, Newman C, Buesching CD, Macdonald DW (2013) A multi- 
metric approach to investigate the effects of weather conditions on the 
demographic of a terrestrial mammal, the European badger (Meles meles). 
PloS one 8: e68116. 

14. Post E, Forchhammer MC, Stenseth NC (1999) Population ecology and the 
North Atlantic Oscillation (NAO). Ecol Bull 47: 117-125. 

15. Blank JL (1992) Phenotypic variation in physiological response to seasonal 
environments. In Tomasi TE, Horton T, editors. Mammalian Energetics: 
Interdisciplinary Views of Metabolism and Reproduction. Ithaca, NY: Corn- 
stock Publishing Associate, pp. 186—212. 

16. Owen M, Wells RL, Black JM (1992) Energy budgets of wintering barnacle 
geese: the effects of declining food resources. Ornis Scandinavica 23: 451—458. 

17. Tinkler E, Montgomery WI, Elwood RW (2009) Foraging ecology, fluctuating 
food availability and energetics of wintering brent geese. Journal of Zoology 278: 
313-323. 

18. Boyce MS (1979) Seasonality and patterns of natural selection for life histories. 
The American Naturalist 114: 569-583. 

19. Kunz TH, Wrazcn JA, Burnett CD (1998) Changes in body mass and fat 
reserves in pre-hibernating little brown bats (Myotis kcifugus). Ecoscience 5: 8-17. 

20. Macdonald DW, Newman C (2002) Population dynamics of badgers (Meles meles) 
in Oxfordshire, UK: numbers, density and cohort life histories, and a possible 
role of climate change in population growth. Journal of Zoology 256: 121-138. 

21. Giroud S, Turbill C, Ruf T (2012) Torpor use and body mass gain during pre- 
hibernation in late-born juvenile garden dormice exposed to food shortage. In: 
Ruf T, Bicber C, Arnold W, Millesi E editors. Living in a seasonal world: 
thermoregulatory and metabolic adaptations. Heidelberg: Springer, pp. 481- 
491. 

22. Harlow HJ (1981) Torpor and other physiological adaptations of the badger 
(Taxidea taxus) to cold environments. Physiol Zool 54: 267-275. 

23. Barclay RMR, Lausen CL, Hollis L (2001) What's hot and what's not: defining 
torpor in free-ranging birds and mammals. Canadian Journal of Zoology 79: 
1885-1890. 

24. Sandell M (1990) The evolution of seasonal delayed implantation. Quarterly 
Review of Biology 65: 23-42. 

25. Yamaguchi N, Dugdale HL, Macdonald DW (2006) Female receptivity, 
embryonic diapause, and superfetation in the European badger (Meles meles): 
Implications for the reproductive tactics of males and females. The Quarterly 
Review of Biology 81: 33-48. 

26. Ferguson SH, Elkie PC (2004) Seasonal movement patterns of woodland caribou 
(Rangifer tarandus caribou). Journal of Zoology 262: 125-134. 

27. Macdonald DW (1976) Food Caching by Red Foxes and Some Other 
Carnivores. Zeitschrift fur Ticrpsychologie 42: 170-185. 

28. Post DM, Snyder MV, Fick EJ, Saunders DK (2006) Caching as a strategy for 
surviving periods of resource scarcity; a comparative study of two species of 
Neoioma. Functional Ecology 20: 717-722. 

29. Easterling DR, Meehl GA, Parmesan C, Changnon SA, Karl TR, et al. (2000) 
Climate Extremes: Observations, Modeling, and Impacts. Science 289: 2068- 
2074. 

30. King D (2005) Climate Change: the science and the policy. Journal of Applied 
Ecology 42: 779-783. 

31. Tuljapurkar S (2010) Environmental variance, population growth and evolution. 
J Anim Ecol 79: 1-3. 

32. Yoda K, Sato K, Niizuma Y, Kurita M, Bost C, et al. (1999). Precise monitoring 
of porpoising behaviour of Adclic penguins determined using acceleration data 
loggers. J Exp Biol 202: 3121-3126. 



33. Shepard ELC, Wilson RP, Quintana F, Laich AG, Liebsch N, et al. (2008) 
Identification of animal movement patterns using tri-axial accclcromctry. 
Endanger Species Res 10: 47-60. 

34. Wilson RP, Shepard ELC, Liebsch N (2008) Prying into the intimate details of 
animal lives: use of a daily diary on animals. Endanger Species Res 4: 123-137. 

35. Ellwood SA, Wilson RP, Addison AC (2007) Technology in conservation: a 
boon but with small print. In: Macdonald DW, Service K editors. Key Topics in 
Conservation Biology. Oxford: Blackwcll Scientific, pp. 105-119. 

36. Hays GC, Forman DW, Harrington LA, Harrington AL, Macdonald DW, et al. 
(2007) Recording the free-living behaviour of small-bodied, shallow-diving 
animals with data loggers. Journal of Animal Ecology 76: 183—190. 

37. Ropcrt-Coudcrt Y, Daunt F, Kato A, Ryan PG, Lewis S, et al. (2009) 
Underwater wingbcats extend depth and duration of plunge dives in northern 
Moms bassanus.J Avian Biol 40: 380-387. 

38. Okuyama J, Kawabata Y, Naito Y, Arai N, Kobayashi M (2009) Monitoring 
beak movements with an acceleration datalogger: a useful technique for assessing 
the feeding and breathing behaviors of sea turtles. Endangered Species Research 
10: 39-45. 

39. Viviant M, Trites AW, Rosen DAS, Moncstiez P, Guinet C (2009) Prey capture 
attempts can be detected in Steller sea lions and other marine predators using 
acceleromcters. Polar Biol 33: 713—719. 

40. Whitney N, Pratt Jr HL, Pratt TC, Carrier JC (2010) Identifying shark mating 
behaviour using three-dimensional acceleration loggers. Endanger Species Res 
10: 71-82. 

41. Markham A, Trigoni N, Ellwood S, Macdonald DW (2010) Revealing the 
hidden lives of underground animals using magneto-inductive tracking. 
Proceedings of the 8th ACM Conference on Embedded Networked Sensor 
Systems 1—42. 

42. Shamoun-Barancs J, Born R, van Loon EE, Ens BJ, Ooserbeck K, et al. (2012) 
From Sensor Data to Animal Behaviour: An Oystercatcher Example. PloS one 
7: c37997. 

43. Gleiss AC, Wilson RP, Shepard ELC (2011) Making overall dynamic body 
acceleration work: on the theory of acceleration as a proxy for energy 
expenditure. Methods in Ecology and Evolution 2: 23—33. 

44. Macdonald DW, Newman C, Buesching CD, Nouvellet PM (2010) Are badgers 
' Under The Weather 1 ? Direct and indirect impacts of climate variation on 
European badger (Meles meles) population dynamics. Global Change Biology 16: 
2913-2922. 

45. Johnson DDP, Jetz W, Macdonald DW (2002) Environmental correlates of 
badger social spacing across Europe. Journal of Biogeography 29: 411-425. 

46. Kruuk H, Parish T (1981) Feeding specialization of the European badger (Meles 
meles) in Scotland. Journal of Animal Ecology 50: 773-788. 

47. Hofcr H (1988) Variation in resource presence, utilization and reproductive 
success within a population of European badgers (Meles meles). Mammal Review 
18: 25-36. 

48. da Silva J, Woodroffe R, Macdonald DW (1993) Habitat, food availability and 
group territoriality in the European badger, Meles meles. Oecologia 95: 558-564. 

49. MuldowneyJ, Curry JP, O'Keeffe J, Schmidt O (2003) Relationships between 
earthworm populations, grassland management and badger densities in County 
Kilkenny, Ireland: The 7th international symposium on earthworm ecology- 
Cardiff- Wales- 2002. Pedobiologia 47: 913-919. 

50. Woodroffe R (1995) Body condition affects implantation date in the European 
badger, Meles meles. Journal of Zoology 236: 183—188. 

51. Newman C, Zhou Y-B, Buesching CD, Kaneko Y, Macdonald DW (2011) 
Contrasting sociality in two widespread, generalist, mustclid genera, Meles and 
Maries. Mammal Study 36: 169-188. 

52. Fowler PA, Raccy PA (1988) Overwintering strategies of the badger, Meles meles, 
at 57-degrees-N. Journal of Zoology 214: 635-651. 

53. Thornton PS (1988) Density and distribution of Badgers in south-west England- 
a predictive model. Mammal Review 18: 11-23. 

54. Kaneko Y, Newman C, Buesching CD, Macdonald DW (2010) Variations in 
Badger [Meles meles) Sett Microclimate: Differential Cub Survival between Main 
and Subsidiary Setts, with Implications for Artificial Sett Construction. 
International Journal of Ecology 859586: 10. 

55. Holmstrup M, Zachariassen KE (1996) Physiology of cold hardiness in 
earthworms. Comparative Biochemistry and Physiology Part A: Physiology 
115: 91-101. 

56. Nuutinen V, Butt KR (2009) Worms from the cold: Lumbricid life stages in 
boreal clay during frost. Soil Biology and Biochemistry 41: 1580-1582. 

57. Macdonald DW (1983) Predation on earthworms by terrestrial carnivores. In: 
SatchcllJE editor. Earthworm Ecology. London: Chapman and Hall. pp. 393- 
414. 

58. Kowalczyk R, Jcdrzcjcwska B, Zalewski A (2003) Annual and circadian activity 
patterns of badgers (Meles meles) in Bialowicza Primeval Forest (eastern Poland) 
compared with other Palaearctic populations. Journal of Biogeography 30: 463- 
472. 

59. Boyle K, Whelan J (1990) Changes in the diet of the badger Meles meles L. from 
autumn to winter. The Irish Naturalists' Journal 23: 199-202. 

60. Revilla E, Palomarcs F (2002) Spatial organization, group living and ecological 
correlates in low-density populations of Eurasian badgers, Meles meles. Journal of 
Animal Ecology 71: 497-512. 

61. Do Linh San E, Ferrari N, Weber JM (2007) Socio-spatial organization of 
Eurasian badgers (Meles meles) in a low-density population of central Europe. 
Canadian Journal of Zoology 85: 973-984. 



PLOS ONE | www.plosone.org 



8 



January 2014 | Volume 9 | Issue 1 | e83156 



Autumnal Activity of the European Badger 



62. Stewart PD, Macdonald DW, Newman C, Cheeseman CL (2001) Boundary 
faeces and matched advertisement in the European badger (Meles meles): a 
potential role in range exclusion. Journal of Zoology 255: 191-198. 

63. Kilshaw K, Newman G, Bucsching CD, Bunyan J, Macdonald DW (2009) 
Coordinated Latrine Use by European Badgers, Meles meles: Potential 
Consequences for Territory Defense. Journal of Mammalogy 90: 1188—1198. 

64. Bucsching CD, Stopka P, Macdonald DW (2003) The social function of allo- 
marking in the European badger [Meles meles). Behaviour 140: 965-980. 

65. Dyo V, Ellwood SA, Macdonald DW, Markham A, Mascolo C, et al. (2010) 
Evolution and sustainability of a wildlife monitoring sensor network. Proceeding's 
of the 8th ACM Conference on Embedded Networked Sensor Systems: 127— 
140. 

66. Dugdale HL, Macdonald DW, Pope LC, Burke T (2007) Polygynandry, extra- 
group paternity and multiple-paternity litters in European badger {Meles meles) 
social groups. Molecular Ecology 16: 5294—5306. 

67. Kowalczyk R, Zalewski A, JPdrzcjewska B (2004) Seasonal and spatial pattern of 
shelter use by badgers Meles meles in Bialowieza Primeval Eorcst (Poland). Acta 
Theriologica 49: 75-92. 

68. Palphramand KL, Newton-Cross G, White PCL (2007) Spatial organization and 
behavior of badgers {Meles meles) in a moderate-density population. Behavioral 
Ecology and Sociobiology 61: 401-413. 

69. Byrne AW, O'Kceffc J, Green S, Slceman DP, Corner LAL, et al. (2012) 
Population Estimation and Trappability of the European Badger (Meles meles): 
Implications for Tuberculosis Management. PloS one 7: c50807. 

70. Markham A, Trigoni N, Macdonald DW, Ellwood S (2012) Underground 
Localization in 3-D Using Magneto-Inductive Tracking. IEEE Sensors Journal 
12: 1809-1816. 

7 1 . Macdonald DW, Newman C, Nouvellet PM, Buesching CD (2009) An analysis 
of Eurasian badger (Meles meles) population dynamics: Implications for regulatory 
mechanisms. Journal of Mammalogy 90: 1392-1403. 

72. Savill P, Pcrrins C, Kirby K, Fisher N (2010) Wytham Woods: Oxford's 
Ecological Laboratory. Oxford: Oxford University Press. 

73. McLaren GW, Thornton PD, Newman G, Bucsching CD, Baker SE, et al. 
(2005) High rectal temperature indicates an increased risk of unexpected 
recovery in anaesthetized badgers. Veterinary Anaesthesia and Analgesia 32: 
48-52. 

74. Thornton PD, Newman G, Johnson PJ, Bucsching CD, Baker SE, et al. (2005) 
Preliminary comparison of four anaesthetic techniques in badgers {Meles meles). 
Veterinary Anaesthesia and Analgesia 32: 40-47. 

75. Kruuk H, Conroy JWH, Moorhouse A (1987) Seasonal reproduction, mortality 
and food of otters (Lutra lutra L.). Shetland Symposia of the Zoological Society 
of London 58: 263-278. 

76. Macdonald DW, Newman C, Stewart PD, Domingo-Roura X, Johnson PJ 
(2002) Density-dependent regulation of body mass and condition in badgers 
{Meles meles) from Wytham woods. Ecology 83: 2056-2061. 



77. Wilson RP, White CR, Quintana F, Halsey LG, Liebsch N, ct al. (2006) Moving 
towards acceleration for estimates of activity-specific metabolic rate in free-living 
animals: the case of the cormorant. J Anim Ecol 75: 1081-1090. 

78. Burnham KP, Anderson DR (2002) Model Selection and Multi-Model 
Inference: A Practical Information-Theoretic Approach. New York: Springer 
Verlag. 

79. Burnham KP, Anderson DR, Huyvaert KP (2011) AIC model selection and 
multimodcl inference in behavioral ecology: some backgrounds, observations, 
and comparisons. Behavioral Ecology Sociobiology 65: 23—25. 

80. Anderson DR (2008) Model based inference in the life sciences: a primer on 
evidence. New York: Springer. 

81. Gbransson G (1983) Denning activity in Swedish badgers, Meles meles. Acta 
Zoologica Fennica 174: 179-181. 

82. Kruuk H (1989) The Social Badger. Oxford: Oxford University Press. 

83. Goszczyhski J, Jedrzejcwska B, Jcdrzejcwski W (2000) Diet composition of 
badgers (Mclcs meles) in a pristine forest and rural habitats of Poland compared 
to other European populations. Journal of Zoology 250: 495—505. 

84. Johnson DDP, Baker S, Morccroft MD, Macdonald DW (2001) Long-term 
resource variation and group size: A large— sample Held test of the Resource 
Dispersion Hypothesis. BMC Ecology 2: 1. 

85. Webb DR, KingJR (1984) Effects of wetting on insulation of bird and mammal 
coats. Journal of Thermal Biology 9: 189-191. 

86. Kimmcl E, Arkin H, Broday D, Berman A (1991) A model of evaporative 
cooling in wetted hide. Journal of Agricultural & Engineering Research 49: 227- 
241. 

87. Kruuk H, Parish T (1985) Food, food availability and weight of badgers (Meles 
meles) in relation to agricultural changes. J Appl Ecol 22: 705—715. 

88. Lima SL (1998) Stress and decision-making under the risk of predation: recent 
developments from behavioral, reproductive, and ecological perspectives. Adv 
Stud Behav 27: 215-290. 

89. Lcnart EA, Bowyer RT, Hocf JV, Ruess RW (2002) Climate change and 
caribou: effects of summer weather on forage. Can J Zoo 80: 664—678. 

90. Clark CW (1994) Antipredator behavior and the asset-protection principle. 
Behavioural Ecology 5: 159-170. 

91. Woodroffc R, Macdonald DW (1995) Fcmalc-fcmalc competition in European 
badgers, Meles mclcs: effects on breeding success. J Anim Ecol 64: 12-20. 

92. Roff DA (2002) Life History Evolution. Sunderland: Sinauer Associates Inc. 

93. McDermott R, Fowler JH, Smirnov O (2008) On the Evolutionary Origin of 
Prospect Theory Preferences. The Journal of Politics 70: 335-350. 

94. Fcro O, Stephens PA, Barta Z, McNamaraJM, Houston Al (2009) Optimal 
Annual Routines: New tools for conservation biology. Ecological Applications 
18: 1563-1577. 

95. Newman C, Macdonald DW (2013) The Implications of climate change for 
terrestrial UK Mammals. Terrestrial biodiversity Climate change impacts report 
card Technical paper. 



PLOS ONE | www.plosone.org 



9 



January 2014 | Volume 9 | Issue 1 | e83156 



